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Q . Abstract 

O . We investigate the baseline range in which the 623 degeneracy in neutrino oscihation probabihties 
(D . 

Q I is absent for fixed values of ^13 and CP violation phase 6cp- We begin by studying sensitivities 

' of neutrino oscillation probabilities to ^13, 623 and 6cp for very- long-baseline neutrino oscillations. 

We show contour graphs of the muon-neutrino survival probability Pii'^ — > i'^) and the appearance 

, probability P(z^e — > u^) on the cos 2^23 — sin 2^13 plane for baseline lengths L = 1000, 5000, 10000, 

o 

' and 12000 km. For each baseline length, it is found that P{i^fj, — > f^) is more sensitive to sin 2^13 at 



O 



Oh! 



energies around its local maximum while it is more sensitive to cos 2^23 at energies around its local 
minimum. On the other hand, the appearance probability P^i^e i^fi) is sensitive to sin 2^13 and 



QhI cos 2^23 only near its local maximum. We observe that the ^23 degeneracy in -P(i^^ i^^) is absent 



at energies around the local maximum of this probability, provided ^13 is sufficiently large. The 
023 degeneracy is also absent in general near the local maximum of P(i^e '^fj,)- Using analytic 



^ I approximations for neutrino oscillation probabilities, we demonstrate that the above observations 

' for L = 1000, 5000, 10000, and 12000 km are in fact valid for all distances. The implications of 

these results on probing the octant of ^23 are discussed in details. 
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I. INTRODUCTION 



The understanding of neutrino masses and mixing matrix is crucial to unveil the mystery 
of lepton flavor structures. The updated SK analysis of the atmospheric neutrino data gives 

Q 

1.5 • 10"^ eV^ < lAm^^l < 3.4 • 10"^ eV^ sin^ 2^23 > 0.92. (1) 

This is a 90% CL. range with the best fit values given by sin^2023 = 1 and Am|]^ = 
2.1 . 10-3 eV^ respectively. An eariie. result based upon LIE analysis gives U 

1.9 ■ 10-3 eV^ < lAm^^l < 3.0 ■ 10"^ eV^ sin^ 2^23 > 0.9. (2) 

at 90% CL. where the best fit values are given by sin^ 2^23 = 1 and Am|i = 2.4 ■ lO'^ eV^ 
respectively. The scenario of za, — »■ z^t oscillation for atmospheric neutrinos has been con- 
firmed by the K2K experiment |3|,|j]. Furthermore the results in the solar neutrino oscillation 
measurements are also confirmed by KamLAND reactor measurements js], ^. Combining 
these measurements, the LMA solution of the solar neutrino problem is established and the 
updated 2a parameter ranges are given by 3] 

7.21 ■ 10"^ eV^ < Am^^ < 8.63 ■ 10"^ eV^ 0.267 < sin^ ^12 < 0.371, (3) 

with the best fit values Am^^ = 7.92 ■ 10^^ eV^ and sin^ 6112 = 0.314. 

Despite the achievements so far in measuring the neutrino mixing parameters, the sign 
of Am|^, the mixing angle ^13 and the CP violating parameter ^cp in the mixing matrix 
remain to be determined. Furthermore, one is keen to resolve the octant degeneracy of 6'23 

i 

The mixing angle ^13 is constrained by the reactor experiments j^, [lo|. The CHOOZ 
experiment |9| gives a more stringent constraint on ^13 with sin^ 26'i3 < 0.1 for a large 
Amg^ (90% C.L.). A recent global fit^based upon three-flavor neutrino oscillation gives the 
2cr upper bound, sin^ 2^13 < 0.124 [7]. It is well known that the mixing angle 61^ can be 
enhanced by the matter effect in Earth. The appearance oscillations — Ve-, ^ii-, 
and the survival mode — * z/^ performed in a very-long baseline have been proposed 
[11] to probe the angle ^13 and the sign of Amg^. Furthermore, the aforementioned very 
long baseline neutrino experiments as well as future atmospheric neutrino experiments are 



proposed to determine the deviation of 6*23 to maximality 



12|. In this work, we focus on 



the mixing angle ^23- We shall provide a global survey on ideal neutrino energies in the 
GeV range and baseline lengths from 10^ km to 10^ km for probing the octant of the mixing 
angle ^23- The muon neutrino survival probability -P(t'^ Vfj) = P^^ and electron neutrino 
appearance probability P(z/e ^n) = PepL are both studied for this purpose. We observe 
that the muon neutrino survival probability P^^ has complementary dependencies on mixing 
angles ^13 and ^23 as the neutrino energy varies. This property is established by studying 
the dependencies of P^^ on cos 26'23 and sin 2^13 while keeping other parameters fixed. The 
choice of the parameter cos 2^23 is appropriate as 

1 1 

-cos2e23 = 2 -sin2^23, (4) 

which is a probe to the deviation of ^23 to the best-fit value tt/A. We find that the de- 
pendencies of P^^ on cos 2^23 and sin 2^13 at energies near local maxima of this probability 
differ drastically from those at energies near local maxima of the same probability. In the 
former case, the probability P^^ is always more sensitive to sin 2^13. Furthermore, the ^23 
degeneracy is absent in this case. In the latter case, the probability P^^ is more sensitive to 
cos2^^23 while the 6'23 degeneracy is generally present. Such information is useful for probing 
the octant of ^23- We also study sensitivities of the probability Pe^ to cos 26*23 and sin 26*13 
with other parameters fixed. We only focus on energies near the local maximum of Pe^ as 
this probability is not sensitive to mixing parameters for energies near its local minimum. 

This paper is organized as follows. In Section II, we compare results on the oscillation 
probability Pe^ obtained by the full calculation with those obtained by various analytic 
approximations. This comparison is essential since analytic approximations will be employed 
for discussions in later sessions. To set up the analytic approximation, we introduce the 
concept of average density which varies with the total neutrino path-len gth inside the Earth. 



Applying full calculations and the two-layer analytic approximations 13|], we identify the 
energy values for local maxima and local minima of neutrino oscillation probabilities P^^ and 
Pe^ for baseline lengths 1000 < L/km < 12000. It is found that the two-layer approximation 
is quite satisfactory compared to the full calculation for computing these energy values. In 
Section III, we first present the dependencies of P^^ and Pe^ on the CP violation phase 
5cp. It will be shown that, unlike Pg^, P^^ is not sensitive to the CP violation phase 5cp- 
We study numerically the effect of CP violation phase to the appearance probability Pg^. 
The result confirms the so-called magic baseline [14, Q, [l^ at L ^ 7600 km where Pe^ is 



rather insensitive to the CP violation phase. After discussions on the CP violation phase, 
we present the contour graphs of probabilities P^^ and Pg^i on cos 26*23 — sin 29i3 plane for 
baseline lengths L = 1000, 5000, 10000, and 12000 km. At all these basehne lengths, we 
shall see that P^^ is more sensitive to sin26^i3 at energies around its local maximum while 
it is more sensitive to cos 26^23 at energies around its local minimum. Such observations 
are then justified by using the two-layer analytic approximations for neutrino oscillation 
probabilities. With this approximation, the baseline lengths and neutrino energies allowing 
an unambiguous determination of 6*23 through measuring P^^ are identified. In Section IV, 
we discuss the prospects of probing the 6*23 octant via measuring Pe^ and P^^^. We then 
conclude in the same section. 



II. THE COMPARISON OF FULL CALCULATIONS AND ANALYTIC APPROX- 
IMATIONS 



We begin the discussions with the relation connecting flavor and mass eigenstates of 
neutrinos, z/^ = Yli Uait^i, with U the Maki-Nakag mixing matrix [17(1 given by 



U 



C12C13 -512013 Si3e 

-S12C23 - Ci2Si3S23e'^CP C12C23 - Si2Si3S23e*^CP 



v 



Cl3'523 

S12S23 - Ci2Si3C23e*'^c;p -C12S23 - Sl2Sl3C23e*'^CP C13C23 



(5) 



where Sij and Cij denote sin^j^ and cos 6ij, respectively. The value for the Dirac type CP- 
phase ^cp ranges from to 27r. The evolutions of neutrino flavor eigenstates are governed 
by the equation 
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where W{t)) = (t'e(i), i^t^it), ^rif)) , Am^^- = mf — rrij is the mass-squared difference between 
the i-th. and j-th mass eigenstates, and V = \/2GFNe is the effegtive potential arising 
from the charged current interaction between v^, and electrons in the medium with N^. the 
electron number density. Numerically V = 7.56 x 10^^^ (p/[g/cm^])ye [eV] with Yg denoting 
the number of electrons per nucleon. We take Ye ~ 0.5 in our calculations. One solves 
Eq. ([6]) by diagonalizing the Hamiltonian on its right hand side. This amounts to writing 
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FIG. 1: The average Earth density along the path traversed by the neutrino as a function of the 
path length L. 

the right hand side of Eq. ([6]) as U' H'U'^i'(t)) with U' the neutrino mixing matrix in the 
matter and H' = diag(i?i, £^2, -Ea) the Hamiltonian after diagonaUzation. To obtain various 
oscillation probabilities described later, we have used the parametrization in for the 
Earth density profile. 

For analytic calculations, we employ the two-layer approximation for the Earth density 
profile Given a path-length L for a neutrino traversing the Earth medium, one can 

divide L into the sum L = Li + L2 + ■ ■ ■ Ln with each Li corresponding to a region with 
a specific matter density. The average density for this path-length is then given by p = 
(piLi + P2L2 + ■ ■ ■pnLn)/L. The Earth medium can be categorized as the Earth mantle 
and the Earth core. If a neutrino only traverses the Earth mantle, we shall use the one- 
density approximation for the analytic calculation with the density defined by the above 
prescription. However, if a neutrino traverses both the Earth mantle and the Earth core, 
one should write the total neutrino path-length as L = 2Lm + Lc with 



L 




(7) 



5 



where R = 6371 km and Vc = 3480 km are the radii of the entire Earth and the Earth core 
respectively while On is the incident Nadir angle of the neutrino. We note that the critical 
Nadir angle for a neutrino to pass the Earth core is 33.17° corresponding to L = 10674 km. 
For L > 10674 km, one separately defines average densities within the path-length and 
the path-length respectively. The average densities as functions of the neutrino path- 
length is shown in Fig. [H For L < 10674 km, there is only one curve for the average density, 
which is represented by the solid line in the figure. Beyond this distance, one can define the 
average density in the core and the average density in the mantle, which are represented 
by dashed and dotted lines respectively. Alternatively, one can also define single average 
density for L > 10674 km by ignoring the distinction between the mantle and the core. This 
is seen from the solid line for L > 10674 km. However, in our analytic calculations, we shall 
adopt the two- density approach for L > 10674 km. 

For analytic calculations, we only compute oscillation probabilities up to the lowest order 
in r = Am2i/Amli. In other words, we set Am2]^=0 in analytic calculations and conse- 
quently the mixing angle Ou and the CP phase 6cp drop out from the oscillation probabilities. 
The probabilities P^^ and Pe/x in the two- layer approximations are given by 19, 20, 21, 22I, 23 1: 

Paj/x = cos^ 6^23 + {u^ + f ^) sin^ 6123 + 2 cos^ ^23 sin^ ^23 {u cost + v sin t) , 
Pe^ = sin2^23(l-w'-y'). (8) 
The quantities u, v and t are defined as 

u = cos(20'")cos(0")-cos(2ej3-2^™)sin(20'")sin(0"), 

V = - cos(2^^^) [sin(0") cos(20'") cos(2^^3 - 2^™) + cos(0") sin(20'")] 



im \ 
I3J' 



where 



+ sin(2^^) sin((/.'=) sin(2^J3 - 29^ 



A m(c) 

im(c) _ 31 Tm{c) 

{Ml,r^^^ = (Ami, + Ar(^) + A-Z'^VS, 

(mL)-(^) = {Ami, + AT^^^ - A™(^V2, (10) 



with 



A^/'^) = J (Amii sin 2^13)2 + {AT'^'^ - Ami, cos 2^13)2. (11) 
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FIG. 2: A comparison of Pe^ obtained by the full numerical calculation and various approximations. 
The thick solid curve denotes the result by the full numerical calculation. The dotted-dashed 
curve denotes the result by setting lS.m'21 = in the full numerical calculation. The dashed 
curve represents the result obtained by the two-layer approximation in the leading order of r = 
Arrigi/Amg]^. The dotted curve denotes the result obtained by one-density approximation in the 
leading order of r while the thin solid curve is that obtained by the one-density approximation in 
the next-to-leading order of r. 

The superscripts m and c denote quantities defined in the Earth mantle and the Earth 
core respectively. For neutrinos traversing only the Earth mantle, one simply sets = 
0, 2Lm = L m. the above equations and recovers well known expressions for P^^ and in 



the one-density approximation 24 1. 



The accuracy of two-layer approximation is shown in Fig. [2] with a comparison of this 
approximation to the full numerical calculation and other approximations. In the calcula- 
tions, we have assumed the normal mass hierarchy and taken sin26'i3 = 0.3, cos 2^3 = 0, 



5cv = 0, Amj^ = 2.4 x IQ-^ eV^, Amj^ = 8.2 x 10'^ eV^, and tan^ ^12 = 0.39 |25|. This 
set of parameters will be adopted for later calculations unless specific mentioning of other 
choices. This set of parameters differ from the most updated best-fit values quoted right 



after Eq. ([3]). However, both set of parameters give undistinguishable results on P^^ and P^ 



MM 
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FIG. 3: Left panel: the energy at the local maximum of P^^, as a function of L. Right panel: 
energies at local maxima and local minima of -P^/^, as functions of L. 

in the energy range concerned here. A comparison made at L = 11400 km has two purposes. 
First of all, it is known that the series expansion in the parameter r is valid for L/ <^ 10'^ 



(km/GeV) 2J, |26|]. Hence analytic calculations performed at this baseline length test the 
marginal region of the condition LjE^ <^ 10^ (km/GeV). Secondly this path-length implies 
that the neutrino traverses both the Earth mantle and the Earth core. Therefore it is also 
a good test to the two-layer approximation. It is seen that the two-layer approximation, 
unlike the one-layer approximation, reproduces well the peak energies of Pg^i, while it gives 
peak probabilities deviating from those obtained from the full calculation by 15% — 20%. 
We also see that the two-layer approximation agrees well with the full calculation in the 
limit Am2]^ = 0. 

For later analysis, we compute energies at local maxima of Pe^ and those at local maxima 
and local minima of P^^ for the baseline range 1000 < L/km < 12000. The results are 
depicted in Fig. [31 We do not study local minima of Pe^ since their values are not sensitive 
to mixing angles Qy^ and Q2z- It is seen that the analytic approximation is satisfactory for 
computing energies at local maxima and local minima of neutrino oscillation probabilities. 
We point out that the energy curves in Fig. Oare calculated with sin 2^13 = 0.3 and cos 26^23 = 
0. It is found that these curves are not sensitive to the values of sin 26*13 and cos 26*23. 
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FIG. 4: The CP phase dependencies of P^^ and for L = 1000 km, 5000 km, 10000 km and 
12000 km respectively. The probabiUty P^^ is described by the thick cm've while Pe^ is described 
by the thin curve. 

III. CONDITIONS FOR THE ABSENCE OF 023 DEGENERACY IN P^^ AND P^^ 
AT DIFFERENT ENERGIES 

A. The dependencies of P^^ and Pg^t on (5cp 

Before concentrating on ^13 and ^23 dependencies of neutrino oscillation probabilities, 
we first study the CP phase dependencies with the full numerical calculations. It is easily 
seen from Fig. H] that P^^ is not sensitive to the CP phase for all distances displayed. On 
the other hand, Pe^ is rather sensitive to the CP phase for L = 1000 km and 5000 km. 
In order to quantify the CP phase dependence of Pg^i, we study peak values of Pe^, which 
occur at energies described by the curve maxl in Fig. [3] for different baseline lengths. This 
peak value for a specific baseline length depends on the CP violation phase 6cp and we 
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baseline length L [km] baseline length L [km] 



FIG. 5: The difference and the ratio of and P™"^ as functions of the baseUne length. 



denote the maximum and the minimum of this value as and respectively. The 

difference and the ratio of these two values as functions of the baseline length L are shown 
in Fig. [5l It is interesting to note that the ratio P^^/ P^^ increases monotonically with the 
baseline length until L = 7600 km. The ratio begins to decrease for a larger baseline but 
remains larger than 90%. In fact, one can see that P^"^^ and -Pg™™ differ by less than 10% 
for L > 6500 km. We point out that 1 — P^'^/ P^^ reachin g to minimum at L = 7600 km 



confirms the so-called magic baseline for the probability Pe^ l4, ll5|, Il6 |. 



B. The dependencies of P^^ and Pg^j on mixing angles ^13 and ^23 

Having studied CP phase dependencies of oscillation probabilities, we now focus on 61^ 
and 6*23 dependencies. In this study we set the CP phase (5cp equal to zero. The results 
are presented in Fig. [61 It is easily seen that the values of P^^ at its local maximum and 
local minimum depend on mixing angles ^13 and ^23 while only the local maximum of Pg^ 
depends on these parameters. This confirms our earlier comments concerning the left panel 
of Fig. [HI We point out that the differences between solid and dotted curves in Fig. [6] reflect 
the effect of sin 2^13; while the differences between solid and dashed curves there reflect the 
effect of cos 26^23. 

We now present contour graphs of P^^ and Peu_ on the cos 26*23 ~ sin 26*13 plane. The 
range for cos 26*23 is chosen such that sin 26*23 > 0.9 [2|, i.e., —0.316 < cos 26*23 < 0.316; while 
sin^ 26*13 is chosen to be less than 0.1, i.e., sin 26*13 < 0.316. The contour graphs of P^^ at 



10 




FIG. 6: The ^13 and 623 dependencies of P^^ and Pe^ for L = 1000 km, 5000 km, 10000 km and 
12000 km respectively. The probabihty is described by the thick cm've while Pe^ is described 
by the thin curve. 

different baseline lengths are presented in Fig. [71 Except for L = 1000 km, we have shown 
contours of P^^ for energies in the vicinity of both local maximum and local minimum of 
this probability. The contour for the local maximum of P^^ at L = 1000 km is not shown 
since P^^ at this energy and baseline length is not sensitive to mixing angles ^13 and 6*23 . For 
L = 5000 km, 10000 km and 12000 km, it is seen that the contours at local maxima of P^^ 
and those at local minima of P^^ behave rather differently. The former are in general more 
parallel to the cos 26'23-axis while the latter are generally more parallel to the sin 26'13-axis. 
We note that the local maximum (max2) of P^^ at L = 12000 km can vary from 0.9 to a 
much smaller value, 0.45, which is a result of significant matter effects. Similarly, due to 
large matter effects, the local minimum (min2) of P^^ at L = 10000 km can vary from 
to a much larger value, 0.4. We also notice that, at this baseline length, the 6*23 degeneracy 
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FIG. 7: The contour graphs of the muon neutrino survival probability P^^ on the cos 2^23 — sin 2^13 
plane. At L = 1000 km, the local minimum of P^^u on the curve mini occurs at = 1.89 GeV. We 
plot the contour graph of by averaging this probability over an 1 GeV energy range centered 
at the above local minimum. At L = 5000 km, the local maximum of on the curve maxl 
occurs at -E = 4.75 GeV while the local minimum of this probability on the curve mini occurs at 
E = 9.24 GeV. We plot the contour graphs of P^^ in the energy range 4.25 < E/GeV < 5.25 for 
the former case and 8.74 < E/GeV < 9.74. The same type of convention applies to L = 10000 km 
and 12000 km. 

is absent for P^^ > 0.1. In general, such a degeneracy is also absent for energies near local 
maxima of Pn^. However, the probabilities are not sensitive to cos 2^23 in those cases. For 
comparisons, we also present contour graphs for the appearance probability Pe^ at different 
basehne lengths. It is clearly seen that Pg^ is only sensitive to sin 2^13 for most cases. The 
sensitivity to cos26'23 only occurs at very long baseline lengths and large values of sin 26*13. 
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P for E=1.39-2.39 GeV (maxl), L=1000 km 



P forE=5.31-6.31 GeV (maxl), L=5000 km 

en ^ " 




FIG. 8: The contour graphs for the oscillation probability Pg^ on the cos 2^23 — sin 2^13 plane. We 
plot contours of Pe^ at energies near the local maximum (maxl) of this probability. 

For example, at L = 10000 km, P^,^ becomes sensitive to cos 2^23 as sin 2^13 approaches 0.3. 
At L = 12000 km, P^.^ becomes sensitive to cos26'23 when sin26'i3 is greater than 0.2. 



C. A global look at the absence of ^23 degeneracy 

In this subsection, we focus on 6*23 dependencies of P^^ and P^^ for general baseline 
lengths. The two-layer analytic approximations for P^^ and P^fj, will be employed for our 
discussions, and observations in the previous subsection shall be justified. It is instructive 
to rewrite Eq. ([8]) in polynomials of cos 26*23: 

f{y,z) = -ay'' + {a + l3)y+{l-(3), 

g{y,z) = -7(1/ -1), (12) 
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FIG. 9: The coefficient a{z) + f3{z) calculated along the energy curve maxl in the left panel of 
Fig. [3l The values of sin20i3 are taken to be 0.1, 0.2 and 0.3 respectively. 

with f{y, z) = P^^, g{y, z) = Pef,, y = cos 26123 and z = sin 26*13. Furthermore, 

a = — - [{u — cost)^ + (f — sint)^] , 

p = ^{l-u^ -v^) + ^[{u- costf + {v- sintf] , 

7 = '-{l-u^-v'). (13) 

We note that the sin 26*13 dependencies of P^^ and Pg^ reside in quantities u, v, cost and 
sint. These quantities also depend on the basehne length L and the neutrino energy E. 
Hence the coefficients a, j3 and 7 also depend on the baseline length L and the neutrino 
energy E. It is interesting to note that a + /3 = 7. Therefore we have 

P,r = a{y^-l), (14) 

using P^e + P^l^l + P^lT = 1 and P^g = Pen with our choice of 6cp = 0. The contour structure 
of Pefj, is straightforward as g{y,z) is only a linear function of y. Hence no 623 degeneracy 
presents in the contour graphs depicted in Fig. [HI Additionally, the sensitivity of Pg^ to 
cos 26*23 is dg{y, z)/dy = —{a{z) + l3{z)). The coefficient a + P evaluated along the energy 
curve maxl in the left panel of Fig. [3]are plotted in Fig.[9]for sin 2^13 = 0.1, 0.2 and 0.3. For 
sin 26*13 = 0.3, a + (3 reaches to the maximal value, 0.5, for L ^ 10500 km. For sin 26*13 = 0.1 
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FIG. 10: The coefficients —a. and a + /3 evaluated along energy curves in tlie riglit panel of Fig. [3l 
The coefficients are calculated with sin 2^13 = 0.2 and 0.3 on the left and right panels respectively. 
The thick curves denote values of —a. while thin curves denote those of a + /3. For solid and dotted 
curves, the thick curves generally dominate over the corresponding thin ones. For dashed and 
dotted-dashed curves, the thin curves generally dominate over the thick ones. 

and 0.2, a + /? rises quickly as the baseline length L surpasses 10674 km. We note that the 
value of Pe/i is proportional to a + /5. Hence a + /5 shown in Fig. [9] is its own maximal value 
for each baseline length L. 

The contour structure of P^^ can be analyzed through the quadratic polynomial /(y, z) 
in y. If —a ^ a + /?, generally there are two solution curves for /(y, z) = "p compatible with 
the ranges of y and z where p is a given value for P^j^. Let us suppose that z = sin 2^13 



is measured in the future 
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28l | with a central value zq. The two solution curves for the 
equation f{y,z) = p then intersect with the straight line z = sin26'i3 = zq at two points 
{yi, Zq) and (?/2, -Zo)- This is actually what we have seen in Fig. [7]for local minima of P^^. If, 
on the other hand, —a ^ a + /5 or even —a ^ a + /?, there exists only one solution curve 
for the equation f{y, z) = p. This is because that the two points zq) and {y2, Zq) can not 
simultaneously satisfy the constraint —0.316 < y < 0.316 since \yi + y2\ = —{a + l3)/a > 1. 
This is actually what we have seen in Fig. [7] for local maxima of P^^. To justify this 
observation, it remains to show that the coefficient —a dominates over a + /5 at energies 
corresponding to local minima of P^^ while the latter dominates over the former at energies 
corresponding to local maxima of P^^. This is clearly demonstrated in Fig. [10] where the 
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coefficients —a and a + /3 are evaluated along energy curves in the right panel of Fig. [31 
We have calculated the coefficients with cos 2^23 = and sin26'i3 = 0.2, 0.3 respectively. 
We remark that other choices for cos 26'23 do not produce noticeable changes on the energy 
curves where —a and a + /9 are evaluated. 

It is easily seen that always dominates over —a when these coefficients are evaluated 
at energies along mcixl or max2 in the right panel of Fig. [31 In such cases, the ^23 degeneracy 
is absent in the solutions of /(y, z) = p. Namely there exists only one solution curve for 
the above equation. Reversely, —a always dominates over a + f3 when these coefficients are 
evaluated at energies along mini. The situation is slightly more complicated when these 
coefficients are evaluated at energies along min2. In this case —a no longer dominates 
over a + (3 for baseline lengths around lO'* km. In fact, with sin 26*13 = 0.3, a + (3 is even 
larger than —a for 9000 < L/km < 10500. This explains the contour structure of P^^ 
at L = 10000 km (see Fig. [7]) where the straight line z = 0.3 only intersects one equal 
probability curve f{y,z = 0.3)=p. The straight line 2; = 0.2 also behaves the same except 
for a very small p. We reiterate that the range for y = cos 2^23 is —0.316 < y < 0.316 due 
to the constraint sin^ 26*23 > 0.9 [2]. Therefore, given z = zq, the equation f{y, zq) = p could 
have only one solution for y if —{a{zQ) + /3{zo))/a{zQ) > 0.632. In other words, such values 
of —{a{zQ) + P{zo))/(y{zo) lead to the absence of 6*23 degeneracy. In fact, the condition for 
the absence of 6*23 degeneracy is even more relaxed. To see this, let us divide our discussions 
according to the true octant of 6^23- 

1. min2, ^23 < 7r/4 

Since a < and a + (3 > 0, the two solutions for y in f{y, zq) = p are both negative for 
1 — l3{zo) — p > while they have opposite signs for 1 — /5(zo) — p < 0. If the true value of 
^23 is less than 7r/4, i.e., the true value of y is positive, then the experimental measurement 
should give 1 — /3(zo) — p < so that a positive solution for y exists. With 1 — i3{zq) — p < 
0, the two solutions for f{y,ZQ) = p have opposite signs and the negative solution has 
a larger absolute value. The negative solution will violate the constraint —0.316 < y ii 
-{a{zo) + /3{zo))/a{zo) > 0.316. For zq = 0.2, -{a{zo) + f3{zo))/a{zo) > 0.316 is valid 
for 8300 < L/km < 10770. For z = 0.3, the above baseline range is extended to 7410 < 
L/km < 10790. 
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2. min2, 6*23 > 7r/4 

With a true value of ^23 greater than 7r/4, i.e., the true value of y less than zero, the 
value of 1 — (3{zq) — p can either be positive or negative. The condition 1 — /9(zo) — P > 
(<) is equivalent to the condition \y\ < {>) {a{zQ) + P{zq)) / {—a{zQ)). For — (^(^o) + 
P{zo))/a{zo) > 0.316, one must have 1 — P{zo) — p > 0. Hence there exist two negative 
solutions for y. In this case, the corresponding solutions for ^23 are both located in the same 
octant. For 0.316 < —{a{zo) + P{zQ))/a{zo) < 0.632, the spurious solution for y may or may 
not violate the constraint y > —0.316. For — (a(zo) + P{zo))/a{zQ) > 0.632, the spurious 
solution for y must violate the constraint y > —0.316, hence the 6*23 degeneracy is surely 
absent. For sin26'i3 = 0.2, the condition — (a(zo) +/5(-2o))/q^(^o) > 0.632 can not be achieved 
along min2. For sin2^^i3 = 0.3, the above condition is satisfied for 8270 < L/km < 10720. 
For —{a{zQ) + P{zQ))/a{zQ < 0.316, 1 — P{zq) — p can either be positive or negative. For 
1 — P{zq) —p > 0, both solutions for y are negative and satisfying the constraint y > —0.316. 
For 1 — i3{zq) — p < 0, both solutions for y satisfy the constraint —0.316 < y < 0.316. 
However their corresponding ^23 angles are situated in different octants. 

Let us summarize the results obtained in this subsection. The coefficient a + (3 dominates 
over —a for energy values along curves mcLxl and max2 for all baseline lengths. Hence the 
6'23 degeneracy is absent along these energy curves for all baseline lengths. The situation 
along the curve mini is just the opposite, the coefficient —a dominates over a + (3 for all 
baseline lengths. Hence the ^23 degeneracy is present for all baseline lengths in this case. 
The issue of ^23 degeneracy becomes more complicated along min2, which we have discussed 
according to the true octant of ^^23- Along the energy curve min2, the non-degeneracy 
baseline range is larger for the 6'23 < 7r/4 case. 



IV. DISCUSSIONS AND CONCLUSIONS 



We have presented the baselines and energies ideal for probing the octant of ^23 through 
neutrino oscillations. The appearance mode —>■ can be studied in a very long baseline 
with the facility of neutrino factory 29] or the more recent proposed (3 beam jsO]. As said, 
the sensitivity of Pe^ to 6^23 is dg{y,z)/dy = —{a{z) + P{z)) where g{y,z) represents Pe/^ in 
the analytic approximation given by Eq. (fT2ll . The maximal value oi a + (3 for each baseline 



17 



length is shown in Fig. [9l At the magic basehne, L = 7600 km, a + (3 = 0.06, 0.21 and 0.38 
for sin 2^13 = 0.1, 0.2 and 0.3 res p ecti vely . For a sufficiently large sin 2^13 and a baseline 



length close to the magic value [ij, ll5|, ll6[], Pg/x is ideal for probing the octant of ^23- 



The probability P^^ is also relevant in neutrino oscillation experiments with neutrino 
factories. The sensitivity of this probability to 623 is determined by the derivative 

r/P 

r = ^ = -2a cos 2^23 + {a + (3) . (15) 

a cos 26'23 

Since a is negative, the sensitivity r is larger for cos2^^23 > 0, i.e., ^23 < tt/A. For a 
measurement performed around a local maximum of P^^, the sensitivity to 6'23 is completely 
determined by the coefficient a + (3, since the coefficient a is generally rather suppressed in 
this case. Along the energy curve denoted by mcixl, a+f3 peaks at 7480 km for sin 2^13 = 0.2 
and it peaks at L = 7350 km for sin26'i3 = 0.3. The values of a + /5 at those peaks are 0.17 
and 0.32 respectively. Along the curve max2, a + f3 peaks around L = 10750 km for both 
sin 2^13 = 0.2 and 0.3 with values 0.43 and 0.48 respectively. 

For a measurement performed around a local minimum of P^^, the sensitivity to ^23 is 
determined by both coefficients —a and a + (3. Along the energy curve denoted by mini, 
the coefficient —a is always close to unity while the coefficient a + (3 is always suppressed 
for all baseline lengths. It is understood that the magnitude of a + /3 determines the size 
of matter effects. Hence the matter effect is small at energies along the curve mini. The 
suppression of a + (3 compared to —a leads to the 6^23 degeneracy as discussed before. The 
behavior of P^^ along the energy curve min2 is more interesting. If the true value of ^23 
is less than 7r/4, the 6'23 degeneracy from the measurement of P^^ is absent in the baseline 
range 8300 < L/km < 10770 for sin 2^13 = 0.2. The above non- degeneracy baseline range 
extends to 7410 < L/km < 10790 for sin 2^13 = 0.3. On the other hand, if the true 6'23 
is greater than 7r/4, the non-degeneracy baseline range does not exist along the energy 
curve min2 for sin 2^13 = 0.2. For sin26'i3 = 0.3, the non-degeneracy baseline range is 
8270 < L/km < 10720. 

The existence of non-degeneracy baseline range along the energy curve min2 has impor- 
tant implications. It can be seen from Fig. [3] that the curve min2 lies in between curves 
maxl and max2. Since the degeneracy of ^^23 is absent on both maxl and max2 for all 
baselines, it is possible that there exists a non- degeneracy region spanned by ranges of the 
baseline length and the neutrino energy. For example, with sin 2^13 = 0.2 and a true value of 
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TABLE I: The baseline range in which the ^23 degeneracy is absent in the probabihty P^^ for all 
energy values between the curves max2 and maxl. The entry corresponding to ^23 > vr/4 and 
sin 2^13 = 0.2 is left blank since, with such set of parameters, there exists no baseline length where 
the condition for the absence of ^23 degeneracy can be satisfied. 



^23 octant 


sin 2% = 0.2 


sin 26*13 = 0.3 


023 < vr/4 
6*23 > 7r/4 


8550 < L/kni < 10680 


7950 < L/km < 10700 




8450 < L/km < 10680 



6*23 less than 7r/4, the 6*23 degeneracy is absent for 8300 < L/km < 10770 for energies along 
curves max2, min2 and maxl. It is of great interest to investigate if the ^23 degeneracy is 
also absent for any neutrino energy larger than the value on max2 and smaller than that on 
maxl. By taking all these energies into account, we find that the 6*23 degeneracy is absent 
for 8550 < L/km < 10680. For a true value of 6'23 greater than 7r/4 and sin 26*13 = 0.3, 
the ^23 degeneracy is absent for 8270 < L/km < 10720 for energies along curves max2, 
min2 and maxl. However, with all energies between curves max2 and maxl considered, 
we find that the ^23 degeneracy is absent for 8450 < L/km < 10680. The non-degeneracy 
baseline range corresponding to different combinations of ^23 and 6*13 values are summarized 
in Table I. 

It is interesting to compare measurements on Pg^ and P^^ since both oscillations appear 
in experiments with neutrino factories [29]. The main issue for comparison is on the deter- 
mination of the true ^23 value under the assumption that both the sign of Amg^ and the 
value of sin 26*13 are known. Let us begin the discussion with a true value of ^23 less than 7r/4 
and sin 2^13 = 0.2. For L < 8550 km, the appearance mode z/g — > is useful for probing the 
octant of 6*23, in particular for L close to the magic value, 7600 km. However, the survival 
mode —>■ is not as useful since the 6*23 degeneracy is absent only at energies near maxl 
and max2. Concerning the sensitivity to 6*23, we note that the differentiation of Pg^t with 
respect to cos 2^23 is —(a + (3). The value of a + (3 increases with L as shown in Fig. [91 It is 
0.21 at L = 7600 km, and 0.26 at L = 8550 km. For 8550 < L/km < 10680, both z/^ 
and z/^ — >• are useful for probing the octant of 6*23. We note that peak positions of Pe^ 
are mostly around 6 GeV. Hence they overlap with the non-degeneracy energy range of P^^. 
From Eq. (fT5|) and our assumption of 6*23 < 7r/4, we find that P^^ is more sensitive to 6*23 
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as compared to Pe/x for the same neutrino energy. For L > 10680 km, Ue ly^, is again the 
only useful mode for probing the octant of 6*23. 

Let us turn to the case where the true value of 623 is greater than 7r/4 and sin 2^13 = 0.2. 
In such a case, the value for Pg^ is enhanced compared to that with y > 0. Furthermore 
Pe^ is always more sensitive to cos2^^23 as compared to P^^ for the same neutrino energy. 
It is clear that the z/g z/^ appearance mode is more useful for probing ^23 regardless the 
baseline length. Although there exists a baseline range where the ^23 degeneracy is absent 
in P^^ for neutrino energies between curves max2 and mcLxl. However this requires a large 
value of sin 26*13, such as sin26'i3 = 0.3. 

It is essential to remark that the above non-degeneracy baseline range is not sensitive to 
the value of Arn"^^, which we have so far taken to be 2.4 ■ 10~^ eV^. Changing the value 
of Amg^ only shifts the probability curves in Fig. H] and Fig. [6] so that positions for local 
maxima and local minima of these probabilities shift accordingly. However, the maximal or 
minimal values of these probabilities remain unchanged. In other words, although the energy 
curves in Fig. [3] are shifted, the coefficients —a and a + (3 plotted in Fig. [ini which combine 
to form Pe^ and P^^ (see Eq. ( |T2|) ). remain the same. The values of these coefficients as 
functions of the baseline length L then determine the non-degeneracy baseline range. 

In conclusion, we have studied the probabilities Pg^ and P^^ for very long baseline neutrino 
oscillations. We focus on sensitivities of these probabilities to mixing angles 6^13, 6'23 and the 
CP violation phase 6cp- Taking 6cp = as an example, we presented contour graphs of 
Pefj, and P^^ in the sin 2^13 — cos 2^23 plane for baseline lengths L = 1000 km, 5000 km, 
10000 km and 12000 km. The energy values chosen for such studies are in the vicinities of 
either local minima or local maxima of neutrino oscillation probabilities. For each baseline 
length, we have found that P^^ is more sensitive to sin 26*13 at energies around its local 
maxima while it is more sensitive to cos 26*23 at energies around its local minima. On the 
other hand, the appearance probability Pg^ is sensitive to sin 26*13 and cos 26*23 only near its 
local maximum. Such findings have been applied to probe the octant of mixing angle 6*23 
assuming that the angle ^13 and the sign of Arn^^ are known. The appearance probability 
Pg^ is non-degenerate in 6^23- The sensitivity of Pg^ to cos 2^23 is studied for baseline lengths 
from 1000 km to 12000 km. We also studied the sensitivity of P^^ to cos 2^23 for the same 
range of baseline length. We have identified the ranges of neutrino energy and baseline 
lengths where the 6*23 degeneracy is absent. We have pointed out that, for a true value of ^23 
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less than 7r/4 and a baseline length between 8000 and 10000 km, the survival mode z/^ 
is equally good as the appearance mode z/g z/^ for probing the octant of 623. 
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